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T "~ -1 \ Abstract 

> . 

t"**- ■ We compare U(l) lattice gauge theory to an effective model of Maxwell and London 

equations. In the effective model there is only one free parameter, the London penetration 
depth A. It turns out that one can get good agreement between both models if one modi- 
fies the usual definition of magnetic monopole currents in £7(1) lattice gauge theory. This 
I/"") ■ comparison also shows that already at small distances fluctuations of the occuring string are 

0^ | important. Further, we investigate the /3-dependence of the penetration depth and determine 

. the suppression of the monopole condensate in flux tubes. 

1 Introduction 

One of the most fundamental phenomena which have to be explained by Quantum Chromo- 
dynamics (QCD) is the confinement of colour charges. Numerical simulations of QCD on a 
space-time lattice demonstrated the confinement phenomenon but could not yet clarify com- 
pletely the mechanism for the formation of the gluonic flux tube in the QCD vacuum. The 
most promising picture for this mechanism is the dual superconductor model |fj]: Dynamically 
generated colour magnetic monopole currents form a solenoid which squeezes the chromoelectric 
flux between quarks into a narrow flux tube. The energy of this flux tube increases linearly with 
its length and confinement is achieved. 

Magnetic monopoles are topological excitations which condense in the confining phase. At 
strong coupling U(l) lattice gauge theory has a confining phase and the monopole density 
is high, whereas it decreases exponentially above the phase transition to the Coulomb phase 
0. Therefore it can be regarded as a prototype of a confining theory, and we will use it 
for studying some features of the confinement mechanism. U(l) as a subgroup of SU(N) is 
also important for investigating confinement in non-Abelian gauge theories, using the Abelian 
projection approach Q: After partially fixing a gauge one can define "Abelian monopoles" which 
could be responsible for confinement. The nature of these monopoles, however, remains obscure 
due to their dependence on the particular gauge that has been chosen. Interpreted in terms 
of 't Hooft-Polyakov-like monopoles they can possibly be regarded as regular extended objects 
which do have physical relevance [Q] . In any case, the identification of monopoles and their role 
in the microscopic mechanism of confinement is still an interesting open problem. 
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There have been many different approaches for describing the interaction of two charges in 
a confining medium. In the context of lattice gauge theory we would like to mention: 



calculation of flux tube profiles by Wilson-loop - plaquette correlators, 

application of London theory of dual superconductivity for the connection of field strength 
and monopole currents, 

• string picture of confinement (as exhibited in the Hamiltonian formulation of lattice gauge 
theories ||). 

The flux tube connecting two (colour) charges can be made visible by calculating observables 
like the energy or action density ||. For studying the confinement mechanism, however, it 
is also important to measure the electric field strength itself, rather than its square. This 
was done for U(i) j7| and also for SU(2) The dual superconductor model assumes that the 
induced monopole currents obey a dual version of the London equation. Together with Ampere's 
law this determines the electric field distribution, which can be compared to the results of 
lattice simulations. In this way one obtains the "London penetration length" of the confinement 
vacuum. In ref. 0] the persistent monopole current distribution in the presence of a charge pair 
is measured, and the results are interpreted as validity of fluxoid quantization known from flux 
tubes in ordinary superconductivity. A careful analysis of lattice results, however, shows that in 
the usual scheme a discrepancy appears between electric field and magnetic current distributions 
compared to the predictions based on an effective model of Maxwell and London equations. In 
this article we propose a modification of the identification scheme for monopoles. Within this 
scheme we get good agreement for electric flux and monopole current distributions between U(l) 
lattice calculations and the effective model. 

Lattice calculations also indicate that the width of the flux tube increases with the distance 
of sources. At least this is true in the case of finite temperature which is unavoidable when using 
smaller lattices. This means that in addition to the linear or term in the potential there is also 
a Coulomb-like term. In the string picture of confinement this behaviour has been explained by 
string fluctuations ||. We investigate such string fluctuations within the dual superconductor 
picture. 



The effective model of coupled Maxwell-London equations has been solved in ref. [10] in 
the continuum for a straight string. For performing accurate quantitative comparison between 
lattice simulation of £7(1) and an effective model based on the dual superconductor picture, it 
is necessary to take care of lattice effects. Therefore, we solve the effective model on the same 
lattice that we use for Monte Carlo calculations. In this way we gain control over anisotropy 
effects as well as finite size effects (using periodic boundary conditions). It will also be possible to 
consider the finite temperature in our UiX) simulations which is due to the finite time extension 
of the lattice. 

In section 2 of this article we discuss our modification of identifying magnetic monopoles and 
determine electric field and magnetic current distributions around a pair of charges. Finally, 
we show the behaviour of the monopole condensate in external electric fields. In section 3 we 
perform a comparison between U(l) lattice gauge theory and the results of the effective model 
and investigate the dependence of the penetration depth A on the coupling. 



2 



2 Electric field strength and monopole currents in U(l) 



We examine the electric field and monopole current distribution in the strong coupling phase 
in the presence of a pair of electric charges. A Euclidean 8 3 x 4 lattice with periodic boundary 
conditions is suitable to compare lattice calculations and effective model with sufficient accuracy 
and to reveal also tiny effects. U(l) gauge degrees of freedom are represented by link variables 
Ufi(x) = e %e ^ x \ where the phases 0^{x) G (— ir, it] correspond to eaA fl (x) in the continuum limit 
(e is the coupling constant and a the lattice spacing). For the simulation we use the standard 
Wilson action 

s = p Yl [l-cosiy(x)], (1) 

where /3 = ^ is the inverse coupling and 6 fiu (x) G (— 7r,7r] is the phase of the plaquette variable 

£V(*) = e^ (x) = U^(x)U u (x + p,)Ul(x + u)Ul(x). (2) 

Static charges can be represented by Polyakov loops L(f) = n£i U^=o[x = (r, ka)]. We are able 
to determine expectation values of physical observables O around a charge pair at a distance d 
by the correlation function 

{L(0)L*(d)O(x)) 

where the angle brackets denote the evaluation of the path integral. Let us recall now the 
definition of electric fields and explain our modification in identifying magnetic currents. These 
definitions turn out to be essential for the validity of a lattice version of the Gauss law and 
ffuxoid quantization. 

In the continuum the phase 6^ of a Wilson loop around an infinitely small square of size o 2 
in the ^-plane defines the field strength F^ by a 2 eF^ u = 9^. On a finite plaquette, however, 
two (nearly) equivalent plaquette angles with values ±7r(=Fe) would correspond to two drastically 
different values of the field strength = ±7r/(a 2 e). There is a large class of functions which 
fulfil the requirements of the correct continuum limit and 27r-periodicity. Therefore we demand 
that the Gauss law connects the field strength with external charges ±e which we put in the 
vacuum. 

A derivation of the Gauss law shows that the appropriate field strength definition depends 
on the action used. Since the path integral is invariant under translations of link phases O^x), 
we perform a variation of the angle O^f, X4) in the expectation value of a pair of Polyakov loops: 

^/D M ^t( f+) r(r_) = a (4) 

For the Wilson action this results in 



E 

1=1 



(L(r + )L*(rU sin # Di ) 



6 



(9 4 Men;) 

where the sum runs over the six space-time plaquettes including the regarded link O^x), and 
the sign accounts for the direction of the link in the plaquette If we identify 

6 1 

a 3 div£(x)= ]T (±)-sin0 Di , (6) 

i= 1 
(«4W6 n i) 
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we get for the Wilson action the natural Euclidean field strength definition 



a 2 eF ilv = sin 6^. (7) 

From the validity of equ. (||) it can be seen that for the Euclidean metric the expectation values 
of the charge density (Tl| as well as electric fields turn out to be imaginary. The corresponding 
observables in Minkowski space are real as they should. 

DeGrand and Toussaint H suggested to identify magnetic monopoles as ends of Dirac strings 



[12]. If the sum over link phases around an elementary plaquette exceeds ±7r one defines that a 
Dirac string is crossing this plaquette. In the case that the net number of Dirac strings through 
the surface of a three-dimensional cube is nonzero a monopole or antimonopole is located in the 
center of this cube. Whereas the location of a Dirac string is gauge dependent and therefore 
cannot have a physical meaning, the monopole number inside a cube is a gauge invariant quantity. 
This is easily seen by expressing the integer valued monopole current through the sum of six 
gauge invariant plaquette phases which form the boundary of a three-dimensional cube: 

Di = 2™%^, (8) 

□ iGcube(^pcr) 

which is equivalent to counting Dirac strings. Each Dirac string carries a quantized magnetic 
flux of 2ir/e, therefore the magnetic current in terms of the electric coupling reads 

J* OT = ^>. 0) 

We suggest to modify this definition and to use instead of (j8|) and @ the dual Maxwell equations 
for the identification of monopole currents. These read on the lattice 

° 2 F °i = Jm £ Vpv> (1°) 

□ i£cube(^pcr) 

with Fat according to our field strength definition (0). Monopole currents defined in this way 
form closed loops on the dual lattice and obey therefore a current conservation law. However, 
it is no longer true that the monopole charge in each cube is an integer multiple of 2-7r/e. Only 
if a separated monopole is investigated in the classical limit (S — * rain) the flux through a large 
surface enclosing it turns out to be quantized due to sin 6 — > 6. 

With the definitions (g) and @ we determine the distribution of the electric field strength E 
and the magnetic currents J m in the confinement phase {(3 = 0.96) for a pair of static charges ±e 
at distances from la to 4a. For our "measurements" we take 20 000 gauge field configurations 
which have been generated by a standard Metropolis algorithm and are separated from each 
other by 200 complete updates. The errorbars shown in the figures indicate the jack knife error, 
which is around 50 % higher than the naive statistical one. As expected, the monopole currents 
behave like a coil, squeezing the electric field into a narrow flux tube. The magnetic monopoles 
form left-handed currents around the electric flux lines. 

In Fig. [I] we show the electric field and the curl of the magnetic current on the QQ-axis 
for a distance d = 3a of charges. The absolute values of both quantities decrease clearly in 
the middle between charges. This is an indication for string fluctuations as will be discussed 
in the next section. A transverse profile of electric field strength and monopole currents, which 
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Figure 1: Longitudinal profile of the electric field and the curl of the magnetic current on the 
QQ-axis for a distance d = 3a of charges: Both quantities give a decreasing signal in the middle 
between charges (zeta = 0). 
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Figure 2: Transverse profile of electric field strength and monopole currents for distances d = la 
and d = 3a, in the symmetry plane between charges (zeta = 0) : For d = la the flux tube is 
concentrated on the link connecting the charges, for d = 3a it spreads out significantly which 
can also be seen by the decreasing monopole currents. 
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Figure 3: The square of spatial (< J^ ni >) and temporal (< J 2 nt >) monopole links in a 
transverse profile of the flux tube for a distance d = 2a of charges. 



characterizes the thickness of the flux tube, is shown in Fig. ^ for distances d = la and d = 3a 
of charges. 

It is also interesting to look at the expectation values of the squared monopole currents in 
the region of the flux tube. The vacuum of the confinement phase without sources is crowded 
with randomly distributed current loops which are "aligned" by an external electric field. Calcu- 
lations on the Euclidean lattice show that both spatial and temporal components of the squared 
monopole current are suppressed in the flux tube, i.e. they are negative compared to their 
vacuum expectation value. One has to take into account, however, that due to the definition 
(|l0| ) the spatial monopole currents are again purely imaginary, like the electric field strength, in 
the Euclidean metric. Therefore, the sign has to be changed when going back into Minkowski 
space. Consequently, the expectation values of the squared components of the monopole current 
behave differently in regions of electric fields, as shown in Fig. ||: The temporal component 
(= the monopole density) is suppressed compared to the vacuum, whereas the square of the 
spatial components is increased by the alignment process. The situation is the same as for the 
calculation of electric and magnetic energy density ||, where the sign has to be changed for the 
square of the electric field strength. The expectation value of E 2 in the flux tube is positive, B 2 
turns out to be negative with respect to the vacuum. The observed decrease in the monopole 
density shows a further physical connection to a real superconductor, where the Cooper pair 
condensate is suppressed in regions of external magnetic fields. 



3 A Classical Effective Model 

For a pair of static point charges p(f) = e[5(f — r+) — 5(f — r_)] the coupled Maxwell-London 
equations in a dual superconductor read 

dwE(r) = p(f), (11) 
rot^(f) = -J^(r), (12) 
E(r) = \ 2 rotf m (r). (13) 
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There is a discrepancy between the Gauss law and the London equation. Therefore, we have 
to introduce a dual Dirac string [13| carrying a quantum of electric flux 

£(r) = ej^dr 5 (3) {r-r), (14) 

where I is an arbitrary line connecting the charges. Instead of the London equation (12) we 
write 

E(r) = A 2 rot f m (f) + £(r). (15) 

The integral form of this equation is nothing else but the fluxoid quantization. 

One can obtain E{r) and J m (r) in terms of the electric charge density p{r) and the dual 
Dirac string £ (f). The analytic solutions of eqs. (|Ti~|), (12) and ( |i~5| ) take the form 



J m (r) = -rotP(r), 
E(r) = P{r) — grad</>(r), 

where 

-r—Pe(^)-T=i — ppr (17) 
47r \r — r\ 

is the Yukawa potential of the sources and 

p «=i^/* -ww (18) 

describes the polarization of the vacuum. In ordinary electrodynamics polarization is given by 
the density of electric dipoles in the medium, in our case it is provided by monopole currents 
around the string. If we assume the string £{r) to be the straight line connecting the charges, 
the total energy of the electric field E(r) can be calculated analytically as a function of the 
distance d. It contains a linearly rising and in addition a Yukawa potential: 



x 



W(d) = ~z I d 3 rE 2 (r) = ^ d — — V a divergent constant. (19) 
2 J 8tta Ait d 

It should be stressed that for a straight string of infinite length (f + — > — oo, r_ — > +oo) formulae 
( |l"6| ) - ( |T8| ) restore the field and current distribution used in ref . ^ . The solution for the electric 
field takes the simple form E(R) = (e/27r \ 2 )Kq(R/ X) in this case. 

As discussed in the introduction we prefer to formulate the effective model of above equations 
on a lattice of the same size as we used for the Monte-Carlo calculations shown in Figs. |l] to ||[ On 
a spatial three-dimensional lattice the electric field resides on the links, whereas the monopole 
currents are located on plaquettes. By constructing the curl of these quantities one changes 
from links to plaquettes and vice versa. Therefore equ. (15) is defined on links, equ. (^) on 



plaquettes. The two coupled equations can easily be solved iteratively for arbitrary paths of the 
string £ . For comparison to Monte-Carlo calculations of U(l) lattice gauge theory we consider 
the fields and currents on the whole lattice and not only in the symmetry plane between the 
charges. Our results give clear evidence that in order to get reliable values for the penetration 
depth A it is necessary to consider fluctuations of the string £: 
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Figure 4: The "penetration length" A, obtained by a fit of the electric field profile in the 
symmetry plane within the effective model considering just one straight string: The strong 
dependence on the distance of charges reflects the broadening of the flux tube. 

• If we try to fit A in the symmetry plane, using only a straight string, it turns out that 
an increasing distance between the charges from d = la to d = 4a would require strongly 
increasing A-values A(d), see Fig. ^j, which in our case range from 0.4 until above 1.0. This 
is of course against physical expectations because the penetration depth A characterizing 
the confinement vacuum should be only a function of the coupling constant. 

• The U{\) monopole current profiles for d = 3a and d = 4a show a minimum in the middle 
between charges (see Fig. [l]). Within the effective model this cannot be explained if the 
string is just a straight line connecting the charges. 

• The ratio of the electric field strength and the curl of monopole current off the QQ-axis 
is far from being constant (see Fig. || for d = 3a), as the validity of the London equation 
would suggest. 

Therefore, we perform a thermal average over some thousands of static string shapes, where 
strings up to a maximal length Z max = 11a are considered. We use the same periodic boundary 
conditions as in U(l) simulations, strings are therefore allowed to reach their final point over 
the period. The temperature of U{\) simulations is 1/A^a, which means that in our simulation 
it is given in lattice units by T = 1/4. Field configurations for a special path of the string have 
a classical energy 

W = \H E lv ( 2 °) 

x,i=l,3 

and contribute according to their Boltzmann weight exp(— W/T). A fit to all Monte-Carlo 
results of E and J m on the whole lattice for the distances d = la, 2a, 3a and 4a gives A = 0.163a 
for j3 = 0.96 with \ 2 ~ 1- As an example for the excellent agreement we show a comparison of 
U (1) and effective model data for the profiles of electric field and monopole current in Fig. ^ (for 
a distance d = 2a between the charges). Because the statistical errors are very small for small 
distances a quantitative comparison is more significant in these cases. In the figure of the current 
profile we have also plotted the results of a measurement according to the DeGrand-Toussaint 
definition. It can be seen that the currents increase by a factor of around 1.5 at the chosen 
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U(l) effective model 




Figure 5: On the links of the three-dimensional lattice we show the ratio of the electric field 
strength and the curl of monopole current in a plane of the QQ-axis for distances of charges 
d = la (left figure) and d = 3a (right figure): For each link the left box with errorbars shows the 
the result of £7(1) simulation, the right one shows the ratio obtained within the effective model 
with A = 0.163a. 

/3-value of 0.96. Even for a single distance d this definition would not allow to fit at the same 
time the exponential falling off and the absolute values of E and J m by the effective model of 
coupled Maxwell-London equations (as long as sin 9^ is taken for measurement of electric field 
strength, which was shown to satisfy the Gauss law). 

The ratio of the electric field strength and the curl of monopole current turns out to be very 
sensitive on string fluctuations. Within the effective model its strongly varying value can only 
be explained by considering fluctuations of the string where the London relation is replaced by 
the fluxoid quantization (as shown in Fig. |^). Far from the charges this ratio should tend to A 2 
(which is difficult to check because of the decreasing signal in Monte Carlo calculations); in the 
vicinity of the charges it takes significantly larger values. This can be interpreted as increasing 
probability that the string runs along the regarded link. If the ratio gets negative the curl of 
the monopole current has changed its sign which indicates that the "string" contribution to 
the expectation value is dominating. For a distance d = 3a this happens not only on the line 
connecting the charges, but also for the transverse component near the symmetry plane (see 
Fig. |5|). The good agreement of this ratio between U(l) and the effective model also shows that 
the classical energy ( |2(i| ) is a good quantity for measuring the total energy for a single string 
shape. Remaining tiny deviations may be due to the fact that we have considered strings only up 
to a certain length. Strings closing periodically are therefore slightly underestimated. Another 
restriction of our effective model is the lack of time dependence of string fluctuations. 

Finally we investigate the dependence of the fitted penetration length A on the inverse 
coupling (5 (see Fig. [j]). For this purpose we perform Monte-Carlo simulations for seven further 
values of (3 from (3 = 0.6 to (3 = 1.0, where the phase transition to the Coulomb phase takes place. 
In the strong coupling regime the effective model works very well, because string fluctuations 
become less important. At j3 — > 1.0 the penetration length increases dramatically, which is in 
agreement with the physical expectation of a divergence of A at the phase transition from the 
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Figure 6: Transverse flux tube profile (d = 2a, zeta = 0.5a) of electric fields and monopole 
currents in U(l) (E,J m ) and in the effective model (£ mo d,4imod) with the fitted parameter 
A = 0.163a. Errorbars have been omitted because they are smaller than the symbols. The 
monopole currents in DeGrand-Toussaint definition are also plotted (J m dt) for comparison. It 
turns out to be larger by a factor of approximately 1.5, independent on the strength of the 
signal. 



"superconducting" to the normal state. Above the phase transition we are not able to describe 
the field and current distribution of U(l) by means of our effective model which shows itself in 
strongly increasing x 2 values. 

4 Conclusion 

In the confined phase of U(l) lattice gauge theory electric field distributions and monopole 
currents around a pair of electric charges can be very well described by Dirac's extension [jl2| 
of Maxwell's electrodynamics supplemented by the dual London equation. The investigated 
properties of the highly non- linear U(l) quantum system are equivalent to an effective model of 
coupled (linear) Maxwell and London equations with the penetration depth A as the only free 




I I I L 

0.2 0.4 0.6 0.i 

beta 



Figure 7: The fitted penetration length A as a function of the inverse coupling (3. 
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parameter. A necessary condition for the agreement of both models for finite lattice constant a 
is the identification of monopole currents in accordance with Maxwell equations. Further, it is 
essential to consider fluctuations of the dual Dirac string in the effective model. Proceeding this 
way the fields for all distances of charges can be described by a universal penetration depth A. If 
the electric field penetrates into the monopole condensate monopole currents in time direction 
are suppressed and consequently the condensate is partly expelled from the region of electric 
fields. This leads to a repulsion between monopole condensate and electric flux lines and to the 
compression of the flux lines in a tube. Finally we would like to mention that the definition of 
monopole currents by the Maxwell equations can be generalized to non-Abelian gauge theories, 
allowing a gauge independent verification of the dual superconductor model in QCD. 
With pleasure we acknowledge fruitful discussions with J. Greensite, A. Ivanov and S. Olejmk. 
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